TEV BLAZARS: STATUS OF OBSERVATIONS 
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Q ■ ABSTRACT The close relation between ground-based TeV observations and satellite borne 

I ' 7-ray measurements has been important for the understanding of blazars. The observations which 

involve the TeV component in blazar studies are reviewed. 
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1. INTRODUCTION 

Six years after the discovery of the first extragalactic 7-ray source at TeV en- 
ergies (Punch et al. 1992) the efforts to understand the non-thermal emission pro- 
cesses have stimulated the field of ground-based 7-ray astronomy. The experimental 
techniques have evolved in a variety of aspects: starting from contemporary obser- 
vations with several instruments covering a substantial part of the multi- wavelength 
spectrum and by measuring the energy spectra and their temporal behavior. 

Extragalactic TeV astronomy has been inspired by the detection of more than 50 
AGNs between 30 MeV - 30 GeV by the EGRET detector (Thompson et al. 1995) 
on-board the Compton Gamma Ray Observatory (CGRO). Nearly all of the EGRET 
AGNs are radio-loud, flat spectrum radio sources in which the radio emission comes 
predominantly from a core region rather than from the outer lobes. The relatively 
high 7-ray luminosity of those sources implies that the emission is relativistically 
beamed (Blandford & Rees 1978; Blandford & Konigl 1979; Maraschi, Ghisellini & 
Celotti 1992) in the direction of the observer. These 7-ray AGNs with their jets 
aligned with the line of sight of the observer are collectively called blazars, but they 
subdivide into two classes, the flat-speclrum radio quasars and BL Lacertae objects 
(review see Padovani et al. 1997). 



The spectral energy distribution (SED) of the 7-ray emission from individual 
blazars clearly indicates a second component in addition to the synchrotron spec- 
trum at lower energies. The second component could be for example due to the 
synchrotron-self-Compton (SSC) mechanism as predicted by Konigl (1981) and 
Marscher & Gear (1985). It has also been realized, that Comptonization of ex- 
ternal radiation could produce 7-rays in blazars (Dermer & Schlickeiser, & Mas- 
tichiadis 1992; Sikora, Begclman & Rees 1994; Blandford & Levinson 1995; Ghis- 
ellini & Madau 1996). Also pion photoproduction by a proton component of the 
jet (Mannheim & Biermann 1992; Mannheim 1993) has been suggested as source 
of 7-rays in blazars. 

The complete coverage of the second component between MeV-multi-TeV en- 
ergies can only be achieved by using satellite-based and ground-based detectors. 
Without either one of them the spectral coverage of the second component would 
be incomplete and the understanding of the emission process at the source signifi- 
cantly less constrained. 

Several ground-based 7-ray observatories (Whipple, HEGRA, CAT) have par- 
ticipated in the study of blazar emission. The Whipple observatory 7-ray telescope 
has started to constrain the high energy end of the second component for one of 
the EGRET blazars: Mrk 421. The detection of Mrk 421 (Punch ct al. 1992; Petry 
et al. 1996) and subsequent monitoring of this source has shown some unexpected 
and extreme properties of the emission process. 



2. VARIABILITY 



Extreme variability on time scales ranging from days (Kerrick et al. 1995) to 
1/2 of one hour (Gaidos et al. 1996) (Figure 1) constrains the size of a presumably 
spherical emission. 

The rapid variability implies a com- 
pact emission region resulting in 
high opacities. To escape the dense 
emission regions the radiation has to 
be substantially Doppler boosted in 
order to avoid absorption through 
7-7 pair production. Assuming the 
TeV emission is correlated with op- 
tical emission a Doppler factor of 
S > 9 has been derived (Gaidos et 
al. 1996). For the variability of the 
uic'hrs) 6 utco^) emission from a spherical blob rcl- 

Figure 1. Light curves of two Mrk 421 flare events ativistic causality requires that the 
of 1996 May 7 (left) and 1996 May 15 (right), size R of the emission region is lim- 
This figure has been adapted from Gaidos et al. ited to R < cT#/(l + z). 
(1996) 
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Gaidos et al. (1996) have shown that R < 10~ 3 — 10 -4 pc. The variability time 
scale places severe constraints on the geometry and the location of the gamma-ray 
emission region. The small emission region might hint that TeV emission comes 
from close to the base of the jet, but this cannot be firmly concluded. 



2. MULTI- WAVELENGTH OBSERVATIONS 

Studying a class of astrophysical objects with a single source provides a limited 
view of the underlying physics and more sources are needed to either derive infor- 
mation about them in a statistical sense and by studying their differences. The 
discovery of Mrk 501 (Quinn et al. 1996; confirmed by Bradbury et al. 1997) 
and 1ES 2344+514 (Catanese et al. 1998a) by the Whipple observatory 7-ray tele- 
scope (E > 250GeV) has emphasized the significance of the high energy coverage 
of blazars: both, Mrk 501 and 1ES 2344+514 are not detected by EGRET. 
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Figure 2. Multi-wavelength observations of Mrk 
501 (Catanese et al. 1998b): (from top to bot- 
tom): gamma-ray, hard X-ray, soft X-ray RXTE 
PCA, soft X-ray RXTE ASM, U-band optical, 
taken during the period 1997 April 2-20. 



Therefore, it is tempting to specu- 
late that some 7-ray blazars emit 
their maximum power at TeV en- 
ergies rather than in the EGRET 
regime. However, simultaneous 
multi-wavelength observations are 
necessary to study possible cor- 
related spectral variability, be- 
cause those sources are highly vari- 
able. Multi- wavelength campaigns 
on Mrk 421 (Buckley et al. 1996) 
showed some evidence for correla- 
tions between X-ray and TeV emis- 
sion. Unusual flaring activity of 
Mrk 501 starting in 1997 February 
(Breslin et al. 1997) lasting until 
1997 August has inspired numerous 
observations including optical, X- 
ray, hard X-ray, MeV-GeV and TeV 
telescopes, some of them were con- 
temporaneous and can be used for 
searching for multi-wavelength cor- 
relations. Figure 2 shows the result 
of one of the multi-wavelength cam- 
paigns (Catanese et al. 1997; Pian 
et al. 1997) revealing possible multi- 
wavelength correlations: 
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Figure 2 shows the fluxes detected at different wavelengths during an episode 
of flaring activity in April 1997. The analysis of RXTE PCA detectors (Catanese 
et al. f998b) has been included in this plot. These new data strengthen the claim 
that TeV emission and X-ray emission in the 2 keV - 25 keV band are correlated. 
Several further conclusion can be derived from the multi-wavelength campaign in 
April 1997: 

(1) The TeV flaring activity coincides with a strong detection of Mrk 50f by the 
OSSE detector aboard the CGRO between 50 - 470 keV. OSSE has detected only 
a few blazars (McNaron-Brown et al. 1995) whereas Mrk 501 shows the strongest 
flux ever detected from a blazar except a high state from 3C 273 (McNaron-Brown 
et al. 1996). X-ray observations of Mrk 501 by BeppoSAX taken in April 1997 
support that the synchrotron power in blazars can peak at hard X-ray energies of 
100 keV (Pian et al. 1997). 

(2) The EGRET detector provides only an upper limit and this seems to be 
consistent with the interpretation that the maximum energy output of Mrk 501 
peaks at TeV energies and not in the GeV regime, already suggestive from the fact 
that the source was discovered at TeV energies and not by EGRET. 

(3) If the hard X-ray flux and the TeV flux are in fact correlated and this were 
true from the optical through the 7-ray regime, one can derive a lower limit on the 
beaming (Doppler) factor of S > 1.5 (Samuelson et al. 1998). 

Summarizing the results from the multi-wavelength campaigns it becomes evi- 
dent that there are correlations between the X-ray regime from 0.5 keV up to hard 
X-rays 100 keV and the TeV emission. Furthermore, it shows that the synchrotron 
emission of Mrk 501 extends up to at least 100 keV which is in contrast to Mrk 421 
which shows a break in the spectrum at 1.6 - 2.2 keV (Takahashi et al. 1996). Also 
the data shows that the amplitude of the flux variations in the X-ray regime are 
smaller than at TeV energies. 

Mrk 501 has also been detected by the HEGRA experiment, the CAT collabo- 
ration and the telescope array (Aharonian et al. 1997; Djannati-Atai et al. 1997; 
Protheroe et al. 1997). Figure 3 shows the daily rates of the CAT experiment 
located in the Pyrenees, the HEGRA experiment on La Palma and the Whipple 
observatory in Arizona. It can be seen that the fluxes measured by the different tele- 
scopes are correlated. Observations by several instruments at different longitudes 
could contribute significantly to study hourly to daily flux variations. To study 
longtcrm variations over weeks and months avoiding gaps of several days Rauben- 
heimer et al. (1997) have demontrated that atmospheric Cherenkov detectors can 
also be operated during the bright moon period. 

2. ENERGY SPECTRA 

The TeV energy spectra of 7-ray blazars reflect both the physics of the 7-ray produc- 
tion mechanism and possibly differential absorption effects at the source or in the 
intergalactic medium. The flaring activity of Mrk 501 has provided excellent statis- 
tics to calculate energy spectra. Spectra have been published by several groups. 
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FIGURE 3. The relative fluxes of Mrk 501 in April 1997 measured the three different ground- 
based 7-ray experiments CAT, HEGRA and Whipple (Protheroe et al. 1997). 



Aharonian et al. (1997) show a spectrum between 1 TeV up to 10 TeV which can 
be lit by a simple power law with a differential spectral index of 2.49 ± 0.11 ± 0.25 
(including the statistical and systematic uncertainties). 

A spectrum derived over a larger energy range by Samuclson et al. (1998) 
(Figure 4) extends from 260 GeV - 12 TeV which is not well described by a simple 
power-law - a simple power law fit yields x 2 test probability of 2.8 x 10~ 7 . Instead the 
spectrum exhibits significant curvature and can be well fit by a parabolic spectrum: 

3(E) = (8.6±0.3±0.7) x lO^- 2 - 22 ± °- 04± °- 05 -(o.47±o.o7)io gl „(E) m -2 s -i TcV -i^ 

Note, that there is not necessarily a contradiction between the two different 
spectra derived since Aharonian et al. (1997) cover a smaller energy range and give 
larger uncertainties in the spectral index. 

Possible causes of curvature at the lower energies could be related to the intrinsic 
spectrum of the source requiring that the spectrum has to be consistent with the 
EGRET upper limits. A small subset of (Samuclson et al. 1998) nearly contem- 
poraneous TeV/GeV observations show that a curved spectrum is consistent with 
the EGRET upper limits whereas a simple power law is not. Causes of curvature 
at higher energies include partial absorption of 7 rays in the intergalactic medium, 
absorption in the jet by pair production on low energy photons, absorption by pho- 
tons near the source, or an end in the primary particle beam energy. Intergalactic 
absorption by infrared photons has been advocated by a number of authors (Gould, 
R. J. & Schreder G. 1967; Stecker & De Jager 1992; Dwek & Slavin 1994; Biller et 
al. 1995) as a potential cause of such absorption. Therefore, the oberved spectrum 
might be a convolution of the blazar spectrum with an external absorption term. 
By considering the extent of this potential modification of the blazar spectrum, 
upper limits to the density of the intergalactic IR field may be determined. Such 
limits have been recently derived by Biller et al. (1998) in a robust manner using 
data from Mrk 421 and Mrk 501. The curved spectrum of Mrk 501 (Samuclson et 
al. 1998) can still accomodate the limits from Biller et al. (1998). 
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FIGURE 4. The energy spectral distribution of Mrk 501 between 260 GcV - 12 TeV, measured 
with the Whipple 10m gamma-ray telescope (Samuclson ct al. 1998). The stars show the 
differential fluxes derived from 15 hours of observations at small (standard) zenith angles, the 
circles show the fluxes derived from large zenith angle observations (details given in Krcnnrich 
ct al. 1999). 



The energy spectrum of Mrk 421 (Krcnnrich et al. 1999) based on short flares 
has been derived with a similar statistical precision as the spectrum of Mrk 501 in 
Samuelson et al. (1998). The resulting spectrum of Mrk 421 is fit reasonably well 
by a simple power law; J(E) = E~ 2 - 54 ± 03 ± 010 photons m" 2 s" 1 TeV^ 1 

This is in contrast to the spectrum of Mrk 501. In figure 5 the spectra of Mrk 
501 and Mrk 421 are shown in comparison. The spectra are different in shape 
and since Mrk 421 and Mrk 501 have almost the same redshift (0.031 and 0.033 
respectively), the difference in their spectra must be intrinsic to the sources and 
not due to intergalactic absorption, assuming the intergalactic infrared background 
is uniform. 

At GcV energies, Mrk 421 is seen by EGRET (Lin ct al. 1992), albeit weakly, 
whereas Mrk 501 is not (Catanese et al. 1997). For the latter it could be argued that 
in a synchrotron-inverse Compton picture it would appear that both the synchrotron 
and the inverse-Compton peak arc shifted to higher energies leaving the EGRET 
GeV energy sensitivity range in the gap between the peaks. As shown in Fig. 5, 
in the energy range 260 GeV - 10 TeV, the spectrum of Mrk 501 is harder at 
lower energies and shows more curvature than Mrk 421. In fact the latter can be 
fit by a straight line (i.e., pure power-law). This is also consistent with the peak 
inverse-Compton power occuring at higher energies for Mrk 501. Therefore, also 
the detailed energy spectra of Mrk 421 and Mrk 501 concur well with the conclusion 
drawn from the multi- wavelength spectrum: in case of Mrk 501 the TeV emission 
is closer to the maximum power peak of the second component than in case of Mrk 
421. Blazar emission can also reveal important information about the underlying 
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source mechanism through spectral variability. Short term flare spectra of Mrk 421 
have been derived in Krennrich et al. (1999). Figure 6 shows the spectra of Mrk 
421 during several flares: a big flare (2 hours: 7.4 Crab) on May 7 1996 (flare I), a 
flare (30 minutes: 2.8 Crab) on May 15 1996 (flare II) and flares from 3 different 
nights (1.5 hours total: 3.3 Crab) in June 1995 (flare III). 
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FIGURE 5. The energy spectral distribution 
of Mrk 501 between 260 GeV-12 TeV (filled 
circles) compared to the spectrum of Mrk 421 
(open stars) (Krennrich et al. 1999) 



FIGURE 6. The energy spectra of Mrk 421 
during different flaring states (Krennrich et 
al. 1999): "Big Flare" (7.4 Crab) on May 
7 1996, flare (3.3 Crab) from May 15 1996 
and a flare (2.8 Crab) at large zenith angles 
from June 20, 29 and 30 1995 



The average fluxes (Crab units) differ by a factor of 2.6 and the signature of 
the flares are different: flare I shows a clear rise, whereas flare II includes the rise 
and the fall and flare III show a constant emission. The energy spectra show the 
same shape and spectral index for all three flares (Krennrich et al. 1999). There is 
no significant change. Changes in the spectral index with the flux have also been 
studied in great detail by (Aharonian et al. 1998) using Mrk 501 observations from 
the HEGRA array of telescopes. Also, no significant variations of the spectral index 
are apparent in the data. 

ACKNOWLEDGMENTS The invitation to the 3rd Integral workshop is greatly appreciated. 



REFERENCES 

Aharonian, F. A., et al. 1997, A & A, 327, L5 



7 



Aharonian, F. A., et al. 1998, astro-ph/9808296 
Biller, S. D., et al. 1998, Phys. Rev. Lett., 80, 2992 
Blandford, R. D. & Konigl, A. 1979, ApJ, 232, 34 

Blandford, R. D. & Rees, M. J. 1978, in Pittsburg Conf. on BL Lac Objects, ed. A. N. Wolfe 

(Pittsburg University Press), 328 
Blandford, R. D. & Levinson, A., 1995, ApJ,441, 79 
Bradbury, S. M., et al. 1997, A & A,320, L5 
Breslin. A. C. et al. 1997, IAU Circ. 6596 
Buckley, J. H., et al. 1996, ApJ, 472, L9 
Catanese, M., et al. 1997, ApJ, 487, L143 
Catanese, M., et al. 1998, ApJ, 501, 616 

Catanese, M., et al. 1998, Conf. on BL Lac Phenomenon, Turku (Finland), in press 
Dwek, E., & Slavin, J. 1994, ApJ., 436, 696 

Djannati-Atai, A. et al. 1997, Proc. of Towards a Major Atmospheric Cherenkov Detector-V, ed. 

by O. C. De Jager, 21 
Dermer, C. D., Schlickeiser, R., & Mastichiadis 1992, A., A & A, 256, L27 
Gaidos, J. A., et al. 1996, Nature, 383, 319 
Ghisellini, G. & Madau, P. 1996, MNRAS, 280, 67 
Gould, R. J. & Schreder G. 1967, Phys. Rev., 155, 1408 
Kerrick, A. D., et al. 1995, ApJ, 438, L59 
Konigl, A., 1981, ApJ, 243, 700 
Krennrich, F., et al. 1997, ApJ, 482, 758 
Krennrich, F., et al. 1999, ApJ, 511, in press 
Lin, Y. C, et al. 1992, ApJ, 390, L49 
Mannheim, K. & Biermann, P. L., A & A, 253, L21 
Mannheim, K. 1993, A & A, 269, 67 

Maraschi, L., Ghisellini, G., & Celotti, A. 1992, ApJ, 397, L5 
Marscher, A., & Gear, W. 1985, ApJ, 198, 114 
McNaron-Brown, K. et al. 1995, ApJ, 451, 575 
McNaron-Brown, K. et al. 1996, ApJ, 474, L85 
Mattox, J. R. et al. 1993, 410, 609 

Padovani, P., 1997, Proc. of XXXIInd Rencontres de Moriond, Les Arcs, Savoie, ed. by Giraud- 

Heraud & Tran Thanh Van, 7 
Petry, D. et al. 1996, A k, A,311, L13 
Pian, E. et al. 1997, ApJ, 492, L17 

Protheroe, R. J. et al. 1997, Proc. of the 25th ICRC (Durban, South Africa), 8, 317 
Punch. M., et al. 1992, Nature, 358, 477 
Quinn, J., et al. 1996, ApJ, 456, L63 
Samuelson, F. W. 1998, ApJ, 501, L17 

Sikora, M., Begelman, M. C, & Rees, M. J. 1994, ApJ, 421, 153 
Stanev, T. & Franceshini, A. 1998, ApJ, 494, L159 

Stecker, F. W., De Jager, O. C. & Salamon, M. H. 1992, ApJ, 390, L49 
Takahashi, T. et al. 1996, ApJ, 470 , L89 
Thompson, D. J. et al. 1995, ApJ Suppl., 101, 259 



8 



